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Abstract 

We study the dynamics of a glassy model with infinite range interactions 
externally driven by an oscillatory force. We find a well-defined transition in 
the (Temperature- Amplitude-Frequency) phase diagram between (i) a 'glassy' 
state characterized by the slow relaxation of one-time quantities, aging in two- 
time quantities and a modification of the equilibrium fluctuation-dissipation 
relation; and (ii) a 'liquid' state with a finite relaxation time. In the glassy 
phase, the degrees of freedom governing the slow relaxation are thermalized 
to an effective temperature. Using Monte-Carlo simulations, we investigate 
the effect of trapping regions in phase space on the driven dynamics. We find 
that it alternates between periods of rapid motion and periods of trapping. 
These results confirm the strong analogies between the slow granular rheology 
and the dynamics of glasses. They also provide a theoretical underpinning to 
earlier attempts to present a thermodynamic description of moderately driven 
granular materials. 

*Present address: LPTHE, 4 Place Jussieu, F-75252 Paris Cedex 05, France. 
LPTENS 00/34, LPTHE 00/36, LPENSL-TH-13/2000. 
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1 Introduction 



In recent years, granular matter has received a growing attention from the physics 
community [p]] . The study of powders is relevant not only because granular materi- 
als have many industrial applications but also because it raises many fundamental 
questions of physical interest. 

We call granular matter all many-body systems constituted by grains of macro- 
scopic size. The grains interact via repulsive dissipative forces due to inelastic colli- 
sions and static friction. The thermal energy scale is totally negligible with respect 
to the typical gravitational energy. In the absence of external perturbations, each 
metastable configuration has an infinite life-time and thermal averaging is meaning- 
less. The static properties of such systems are hence very interesting, the sandpile 
problem being the paradigm |]]. 

Powders flow only when energy is supplied externally. This can be done by 
applying a shear or a vibration. The dynamics of granular matter presents a very 
rich phenomenology that depends not only on the intensity of the drive, generically 
called T, but also on the way the granular system is driven A weak driving 
force can be provided by applying 'taps' to the systems, as has been done in the 
pioneering experiments of the Chicago group |J, where the parameter T is the 
reduced acceleration of the taps. Recently, Nicolas et al. investigated the dynamics 
of a powder by sinusoidally shearing it in a weak manner j3|; V is here the maximal 
amplitude of the strain. Experiments reveal that the time evolution, in the gently 
driven situation, is characterized by an extremely slow dynamics |2|, |3], |J. In the 
tapping experiments for instance, the density still evolves after 10 5 taps 0. When 
the energy injection is much stronger, the granular matter eventually becomes fluid, 
it behaves essentially like a dissipative gas and it is described by a hydrodynamic 
theory that takes into account energy dissipation through the collisions between the 
grains [0], ||. In this work we shall focus on the gently driven regime and we shall 
not address the strongly perturbed situation. 

This phenomenology is clearly reminiscent of the behavior of glass forming sys- 
tems for which the control parameter is the temperature T (or the density p) 
@ 0) H> El> 113 • At high T (say), the samples are in the liquid or gaseous phases. 



When T decreases, their dynamics becomes exceedingly slow, and may even appear 
completely stopped during the experimental time-window. However, at temper- 
atures above but close to the glass transition the relaxation reaches a stationary 
regime characterized by the decay of all correlations in two steps, the second decay 
being related to the structural relaxation. At temperatures below the glass tran- 
sition temperature, the structural relaxation time t r depends on the time spent in 



the glassy phase (the 'waiting time' t w , or 'age') and typically t r oc t w [|TTJ. A 
stationary regime cannot be reached experimentally. This is the aging effect which 
has been observed in a wide spectrum of glassy systems such as plastics [O], spin- 



glasses |13| , glycerol |14j , dielectric glasses [15j , complex fluids [0] , phase separating 



systems |I7|], etc. 

Recently, the similitude in the dynamics of granular matter under vibration and 
glass forming materials has been rationalized by Liu and Nagel [?[]. These authors 
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proposed a phase diagram that unifies the physics of glassy systems and granular 
materials. In its simplest version, the diagram has three axis (T, p, T). The (T, 
p) plane describes the physics of glasses, while the (p, Y) plane describes the one 
of athermal driven systems, like powders or foams. In the low-T, high-p, small-r 
region the system is generically jammed or presents glassy features. In this work, 
we focus on the (T, T) plane of this phase diagram. The drive axis V can represent 
two types of forces: (1) 'shear-like' forces that do not derive from a potential and 
hence do work on the sample (2) 'tapping-like' forces that do derive from a 
potential but do work on the samples when they depend on time. Both modify the 
dynamic behavior and the goal is to identify how, and to which extent, in a general 
manner. One possible scenario is that the age of the driven system stabilizes at 
a power dependent level, typically t T oc r _1 . In the rheological language, this is a 
shear-thinning behavior. Some examples are given by domain growth under flow [pl|] 
or by the non-linear rheology of complex fluids |[20|| , and it has been captured by a 
number of models |[T, |23[. In particular, a (Shear/Temperature) phase diagram 
for glassy systems has been derived in Ref. plR. Another scenario, realized in the 



present paper, is that aging is not stopped, at least in a well-delimitated region of 
the phase diagram. 

The behavior of moderately driven granular matter has received a lot of exper- 
0, |, [H |5f, numerical % ||, |7|, g|, g| and theoretical @, 0, 0, 



imental 



lT| |37], attention. All these studies have demonstrated the glassy 
nature of granular compaction below a critical amplitude of the drive, T*. This has 
been first revealed by the very slow relaxation of the density, but memory experi- 
ments |24]] and simulations ^ZE, |31| inspired by earlier spin glass studies (11], [13| have 
also given support to this conclusion. 

The relation between granular matter and glassy systems is widely assumed. 
Indeed, many of the models which have been proposed to study granular compaction 
are directly adapted from glassy models || |26|, |27|, [28], [!7|, Q . Usually, the 

drive T in granular matter is related to the temperature T in glasses. However, 
the assumption that T = T(T) is a highly non-trivial statement and there is, to 
our knowledge, no microscopic approach that justifies it. We do not make such an 
assumption here. In this respect, Mehta et al. M have built a phenomenological 
two-temperature stochastic model based on the observation that the slow granular 
compaction is basically a two-step process: in this model, a short-time process stands 
for the fast independent-particle relaxation, while a slow one stands for cooperative 
rearrangements. The recent experiment of Nicolas et al. || clearly proved the 
existence of these two (uncorrelated) processes. A two-step process, each thermalized 
with its own temperature is precisely the output of previous studies of the constantly 
driven dynamics of glassy systems in a thermal bath (TBI, BT 
numerical simulation of a sheared supercooled liquid |39[ . 



confirmed by the 



In this paper we study the dynamics of a glassy system permanently perturbed 
by a time-dependent Hamiltonian force P0| . Our aim is to identify which properties 
correspond to those observed experimentally in granular systems and whether an 
effective temperature for the slow degrees of freedom is generated in this weakly 
athermal system. In some sense, our approach is 'orthogonal' to previous ones. 
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We do not propose a new model for the slow granular rheology, but rather ask the 
following simpler question: What is the behavior of a glassy system subjected to 
periodic driving forces? To answer this question, the (T, T, to) phase diagram is 
explored. This is done by studying the dynamics of the, by now standard, mean- 
field glass model, the p-spin glass model 0, £43], f44|, under a time-dependent driving 
force. [Albeit its name, that follows from historical reasons, this model represents a 
structural glass and not a spin glass.] We then discuss how our results have to be 
interpreted in the context of granular materials. 

The paper is organized as follows. In the next Section, we recall the definition 
and main properties of the p-spin glass model. Its behavior under an oscillatory 
drive, for the choice of parameters p = 2 and p > 3, is examined in Sections |3] 
and ^, respectively, in two ways: analytical and numerical solutions of the spherical 
model, and Monte-Carlo simulations of the Ising model. We discuss our results in 
the context of slow granular rheology in Section |5|. Our conclusions close the paper 
in Section ||. 



2 Definition and zero-drive behavior of the model 

The spherical p-spin glass, when p = 2, is simply the spherical version of the 
Sherrington-Kirkpatrick spin glass and it was introduced by Kosterlitz et al. JRj] 
as an exactly solvable model. It is in fact equivalent to the O(n), n — > oo, model 
for ferromagnetic domain growth in three dimensions. Its statics [£|f| and its relax- 
ational dynamics have been extensively studied [46], [47], [RJ. The p > 3 spherical 
models instead are simplified models for supercooled liquids and glasses, in the sense 
that they give a theoretical framework to understand both the statics [11, 43, [4£j] 
and the non-equilibrium dynamics |50| of glass forming systems. 



We study these two cases as generic glassy models, taking advantage of the fact 
that exact equations for their driven dynamics can be derived. 



2.1 Model 

The p-spin glass model in its spherical |^4| and Ising [[E| [51|] versions is defined by 
the Hamiltonian 



-fO[ s ] — ^ Jiii 2 ...i p Si 1 Si 2 . . . S{ p , (l y 



«l<«2<-<i P 

where the couplings J% x i^...i p are random Gaussian variables with zero mean and 
variance (Ji^...^) 2 = p!/(2A^ p_1 ). The spins {si,i = 1, • • • , N} may satisfy a global 
spherical constraint J2i s i 2 (t) = N, or be Ising variables Si = ±1. In the spherical 
version, the driven dynamics is given by the following Langevin equation 

^r = -sS)-^ s -w +/ ™ , ™°w + ''-w- < 2 > 
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where the parameter z(t) ensures the spherical constraint and rji(t) is the thermal 
noise, taken from a Gaussian distribution with zero mean and variance Ik^T . [In 
what follows we use units such that the Boltzmann constant k-Q is one.] The thermal 
bath temperature, T, can possibly be zero. In the Ising version one can still write 
down a continuous Langevin equation by using soft spins and taking the Ising limit 
at the end of the calculations. The spherical version is simpler to treat analytically 
while the Ising version is much simpler to deal with using Monte-Carlo simulations, 
since the spins are bimodal variables. 

In order to mimic tapping experiments, a periodic time-dependent Hamiltonian 
force fj AJ>vmG (t) has to be added to the right hand side of the Langevin equation (g). 
The simplest periodic time-dependence one can think of is a cosine form of period 
r = Itx juj. Thus, we are lead to add a magnetic field in the Hamiltonian (Tip, 



A' 



Hj[s] -> Hj[s] + cos(cut) h Mt) • (3) 

i=i 

In the most realistic numerical experiments, the tapping is modelled by a two-step 
dilation-relaxation process while in more schematic models, the driving force is 
not explicitly time-dependent [f2(| [JT|. Two types of spatial dependence of 

the field will be considered below: constant, hi = h, for all sites i, and random 
hihj = h 2 5ij. There are however no physical differences between the two situations, 
since a spatially constant field is as decorrelated from the ground states of the 
Hamiltonian Hj[s] as a random one. 

2.2 Dynamical equations 

The dynamics of the spherical version of the model is better analyzed in terms of 
the autocorrelation function C(t,t') = (s i (t)s^(i / ))/A^ and the linear response 
function R(t,t') = (si(t)i]i(t')) / (2T N) since, in the thermodynamic limit N — > 



oo, C and R verify closed Schwinger- Dyson equations which read, for t > t', |50 



dC{t,t') = _ z ^ c ^ + P(P 1) f dt"C p - 2 {t,t")R{t,t")C{t",t') 



Of 



+ 2 



P ■* 



dR(t, t' 



[ dt"C p - l (t,t")R(t\t") + h 2 cos(tut) [ dt"cos(cot")R(t',t") 
Jo Jo 

z{t)R{t,t') + p ( p ~^ [ dt"C p - 2 (t,t")R{t,t")R{t",t') , 
2 Jt' 

z (t) = T + — f dt"C p -\t, t")R(t, t") + h 2 cos(cjt) f dt" cos(ut")R(t, t") . 
2 Jo Jo 

(4) 

These integro-differential equations are complemented by the equal-times conditions 
C(t,t) = 1, R(t + ,t) = 1, the symmetry of the correlation, C(t,t') = C(t',t), and 
causality, R(t',t) = 0. In deriving these equations, a random initial condition at 
time t — has been used, which can be interpreted as an equilibrium configuration 
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at infinite temperature. An infinitely fast quench towards the final temperature T 
is performed at t — and the evolution continues at subsequent times in isothermal 
conditions. The energy density e(t) = N~ 1 (Hj[s]) is related to the constraint z(t) 
through 

e(t) = ~ [T - z{t)\ . (5) 

In Eqs. (|), the oscillatory field has been chosen to be constant in space. 

In the p = 2 case, the dynamics simplifies considerably since it can be solved 
directly from the Langevin equation (0). Indeed, this set of N differential equations 
is diagonalized by using the basis of eigenvectors of the symmetric matrix J^. De- 
noting /x the eigenvector associated to the eigenvalue fi, and = /xs and = fi-h 
the projections of the spin and field onto the eigenvectors, one obtains 

ds (t) 

— = {n- z(t))s^t) + cos(c^) + Vlt (t) . (6) 

Here we have considered, for convenience, a spatially uncorrelated random field. 
The autocorrelation and response functions become 

C(t,t')= j dpp^is^s^t')) , 

(7) 



where p(p) = \/4 — p 2 /2ir for p G [—2,2], and zero otherwise, is the density of 
eigenvalues of the random matrix J^-. The spherical condition reads C(t,t) = 1 
and after projection, the amplitude of the field is random with zero mean and 
variance h^h^ = h 2 ^^. 

2.3 Fluctuation-dissipation relation and effective tempera- 
tures 

We shall be interested in the fluctuation-dissipation theorem (fdt) and its possible 
modifications. For driven systems, we do not expect this relation to be satisfied. 
In equilibrium, any correlation function C(t,t') = (0(t)0(t')) [we assume, without 
lose of generality, that (0(t)) = 0], and its associated linear response R(t,t') = 
5(0(t))/ 8f(t')\f=o, where the perturbation / modifies the Hamiltonian of the system 
according to H — > H — fO, satisfy the fdt 

„, A 1 dC(t,t') 

R(t,t') = t>t'. 8 

T dt' ~ y 1 



For non-equilibrium systems, a possible extension is [50, 52 



i\ 1 dC(t,t') 

R(t, t') = ——. — ^—L t>t' . 9 

1 ; T EFF (t,f) dt' ~ 1 ; 



6 



Naively, this equation is simply a definition of the two-time function T EFF (t, t'). This 
extension becomes non-trivial when one realizes that, in the long waiting-time limit 
of many solvable models, this function only depends on times via the correlation 
function itself, 

T EFF (t,t') = T FFF [C(t,t')]. (10) 

By extension, it has been proposed that this exact result for solvable, mean-field like 
models, will also apply to more realistic models with, e.g., finite range interactions. 
Graphically, a convenient way of checking this ansatz is to represent the integrated 
response function £') = dt"R(t,t") as a function of C(t,t') at fixed t' and 
parametrized by the time difference t — t' |53|] . Equation ([TOD implies that a master 
curve (i.e. independent of if) x(C) is asymptotically reached. 



In the aging case, three families of models have been found JTTJ: (1) glassy 
models, for which the x{C) curve is a broken straight line with a first piece, from 
C — 1 to C — (Jea = lim^oo linty^oo C(t, t'), of slope — 1/T and a second piece, 
from C = Qea to C = 0, of slope — 1/T EFF , T EFF being finite; (2) 'domain growth 
models', for which the x(C) curve is still a broken straight line with the second piece 
having T EFF = oo; (3) spin-glass models, for which the x(C) curve has a straight line 
piece from C — 1 to C = (Jea and a curved piece from C = (Jea to C = 0. 

The modification of these plots in a system driven by shear-like forces has been 
studied at the mean- field level in Refs. |1| 21f| . Numerically, the same behavior has 



been found in a sheared supercooled liquid |39fl . Finally, Langer and Liu studied a 
sheared foam and studied the fdt by comparing the fluctuations of the stress on 
the boundary and the corresponding compression modulus 



At the theoretical level, it was shown in Ref. [52] that the factor T EFF is indeed 



a bonafide temperature: it is commonly called 'effective temperature'. A very ap- 



pealing connection of this factor with the ideas of Edwards [38], in the context of 



granular matter, and Stillinger- Weber |55[ in the context of glass forming liquids are 



currently being explored EJ, 57, 58|. We shall come back to this point below. 

Finally we wish to recall a rigorous bound that controls the maximum deviation 
from fdt that can be observed in an out of equilibrium system with Langevin 
dynamics ]59[. When a time- dependent force of period r is applied on the system 



the bound takes the form 
"'"- T (l)C(l.s) 



r r ds a 

Jt w V 



ds 




-TR(t,s)\ < J AM dsD 2 (t,s) — (11 



where W is the total work done by the external force h on the system, per period, 



W = J t T ds (v(s) ■ h(s)), Af is a system dependent numerical factor and D(t, s) is a 
two-time correlation that in most cases of interest is again bounded by a numerical 
constant. In particular for the model we treat in this paper, both M and D 2 (t, s) 
equal one. Hence, for this model we have 



rtw+r / 

/ ds 

Jt w V 



ds 




f<)( ' {LS) -TR{t,s))< dr(-). (12) 
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2.4 Zero-drive dynamics: Jammed states and Reynolds di- 
latancy 

In the absence of a driving force, the model has a dynamical transition at a (p- 
dependent) critical temperature T c . For instance, T c = 1 for p = 2 and T c ~ 0.6123 
for p = 3. Above T c , the equilibration time is finite, and the system reaches 
equilibrium. In this case, both time-translation invariance (tti) and fdt are sat- 
isfied and Eqs. (§) reduce to the mode-coupling equation for the so-called F p _i 
model |TT|, |43], ^(J. The relaxation time diverges at T c , and below T c , the system 
does not equilibrate with its environment. A quench from the high temperature 
phase towards the low temperature phase is followed by the aging dynamics de- 
scribed in the introduction: the relaxation time increases with t w , which means that 
tti is lost. The fdt is modified in the way described in the previous section. One- 
time quantities have a slow relaxation — typically power laws and, for instance, the 
energy density converges to a value e(t — > oo) = e t h{T) which is higher than the 
equilibrium value e eq (T). 



From the static point of view [El, |6T|, T c is the temperature below which the 



energy landscape becomes fractured into many metastable states. 'Many' means 
that their number A/"(T, e) is such that the thermodynamic limit 

exists and is finite. The so-called complexity (or 'configurational entropy') S(T, e) = 
lnjV(T, e) is hence extensive in a range e G [e eq (T), e t h(T)]. The effective tempera- 
ture T EFF defined in the previous Section from the fdt violations can be computed 
directly from the dynamical equations. A remarkable result is that it is also given 
byH 



1 <9E(T, e) 



2~rff de 



(14) 



e=e*h 



at T = 0. [At finite temperature the free energy density replaces the energy density 
in Eqs. (0) and (fT^D-l This relation is very similar to the thermodynamic definition 
of temperature and to the definition of Edwards' compactivity in granular systems. 
We shall come back to these similarities in Section |^. 

The complexity is a useful quantity to understand, for example, the difference 
between the dynamics of the p = 2 and p > 3 models. For p = 2, one has e eq = e t h, 
and the complexity is irrelevant, S = 0. It follows that T EFF = oo. This case falls 
in the 'domain growth family' defined in the previous subsection. Instead, if p > 3, 
e e(? < eth, S > in the region [e eg , e t h], and T EFF is finite. 

If, instead of a quench from a high temperature at t = 0, the dynamics starts 
from one of the metastable states with energy density e eq < e < e t h, the system 
may escape this state, but in a time scale which diverges with N. Ergodicity is 
broken, and the dynamics is blocked. We could call these metastable states the 
'jammed states' of the system. (They are known as the tap states in the spin glass 
literature plfl.) 



S 



When a driving force is applied to the system, the dynamics also depends on 
the initial state. The driven dynamics following a quench under shear-like forces 
has been investigated in Refs. fl8|, [21]]. When starting from a low- lying state, and 
for a small enough driving force, the system remains trapped (jammed). Only by 
applying a large force will the system escape the deep state and have a dynamics 
like the one following a quench [18|. Hence, the system becomes able to move only 
by raising its energy density: this is, in the glassy context, the Reynolds dilatancy 
effect [0 0. 

The effect of these trapping states on the dynamics following a quench are beyond 
the mean- field approximation, Eqs. (|J). We propose to investigate this interesting 
aspect by simulating in Section [4.3| the Ising version of the model keeping finite the 
number N of interacting spins, so that jammed states will be dynamically accessible 
by the system. 



3 Driven spherical Sherrington-Kirkpatrick model 

In this Section we study analytically the driven dynamics of the p = 2 case. 



3.1 Existence of a dynamic transition 

Let us first study the p = 2 case at T = 0. The Langevin equation @ can be 
integrated out; this yields 



l + / dt'cos(ujt')e'^n(t') 
Jo 



where the function 



Q(t) = exp 



dt'z(t' 



o 



(15) 



(16) 



has been defined. The spherical constraint determines Q(t) through the equation 



n(t) 2 = f(t) + h 2 dt' dt"f (t 

with f(t) given by 



V + 1" 



tt(t')tt(t") cos{ut') cos(cut") , (17) 



f(t) = J\p(ri<** 



The autocorrelation and response functions, together with the energy density, read 



C(t, t v 



f 



t + U 



t rt 



+ h 2 dt' dt"cos(cut')cos(cut")n(t')ft(t")f 
Jo Jo 



e(t) 



z{t) _ld\n[Q(t)} 
~ ~2 di 



(19) 
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100 



h = 2.5,a; = 0.6. 




80 100 120 



Figure 1: The function Q(t) exp(— 2t) in the three cases h = and h = 2.5 with 
uj = 0.8 (glassy phase) and 0.6 (liquid phase), from bottom to top. The points are 
the direct integration of Eq. (|TTD, while the full lines are fits to Eqs. (|2"0"|). Except 
at very short times, the approximate forms in Eqs. (|20|) yield the correct behavior. 



Equation (|ITD determines and, consequently, all dynamic quantities. Unfortu- 
nately, it cannot be completely solved analytically except in some asymptotic limits. 

Two special cases have been studied previously, namely the zero field behavior 
of the model, and the case of a dc field. In the former, as discussed in Section [|, 
there is no finite equilibration time below T c — 1 and the system ages forever [[R], 
In the presence of a dc magnetic field of amplitude h [f0|, Eq. (|P7|) yields 



Q{t) ~ exp(A°(/i)t), where X°(h) = (2 + h 2 )/y/l + h 2 . The main effect of the 
field is to introduce a new time scale t r °(h) = [A°(/i) — 2] _1 in the problem. For 
times t w <C t r °(h) after the quench, aging is observed as in the zero-field case, 
whereas at later times t w 3> t r °(h), the system has reached its equilibrium state in 
a field: aging is interrupted. The energy density e(t) converges to eoo(cj = 0, h) — 
-(2 + / i 2 )/(2 v / TT7?). 

At finite o>, the numerical solution of Eq. displayed in Fig. [I] suggests the 
following asymptotic behavior: 

h<h*(u) Q{t) ~ c — (6cos(2o;t + 0) + l) , 

h > h*(uj) n(t) ~ d exp [X(lu, h)t + acos(2cut + 0)] , 

where c, c', a, b < 1 and are numerical constants, and with X(u, h) > 2. This 
means that if the amplitude of the driving force is sufficiently small, the ac field does 
not introduce a new time scale and a full aging behavior is observed. For stronger 
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amplitudes of the field, a time scale defined by 

t r (u,h) = [X(u,h) - 2]- 1 (21) 

is generated and the relaxation time becomes finite: aging is stopped by the driving 
force. There is hence a critical field h* separating these two different regimes. It is 
clear from Fig. [1] that the agreement between Eqs. fl2"0|) and the numerical solution 
is very good, after a short transient. 

In the following we characterize more precisely these two phases as well as the 
transition between them. 



3.2 The transition line in the plane (u),h) 

We have seen that, at fixed pulsation to, there exists a well-defined transition line 
where the relaxation time in an ac field diverges, allowing to distinguish between 
a 'glassy' and a 'liquid' state. We anticipate the discussion of the last section to 
emphasize that this transition is a non-equilibrium phase transition which is hence 
of a different nature that the transition taking place at T c in the absence of the 
driving force. 

The transition line h*(u) may be understood and estimated from a simple phys- 
ical argument. At u = 0, the dc field introduces a finite relaxation time t r (h). The 
most naive requirement for the system to keep a finite relaxation time in an ac field 
of period r is given by t r °(h) < r. The transition line is then estimated by the 
relation 



t r \h* 



UJ oc 



2 + h* 2 
y/l + h* 



(22) 



A more refined computation can also be performed. Using the fact that, when 
h > h*, Q(t) is given asymptotically by Eq. ( p0|) with A > 2 allows us to neglect the 
term f(t) in the right hand side of Eq. (|TTD; thus 



tt(t) 2 ~ h 2 / dfip(ti)e 



3 2 M i 



dt' cos{ut')Sl{t')e 



(23) 



Inserting the form Q(t) = exp(At + B(t)), and formally integrating by parts yields 

2 



.Bit) 



dfip(p) 



3? 



Aujt 



E 

fc=0 



(X(u, h) — fi + iuo) 



k+l 



d k e B ® 
dt k 



(24) 



Since the function B(t) is periodic with angular velocity 2u, this corresponds to 
a development in powers of uj. In the small frequency limit, the equation can be 
closed by keeping only the leading terms of the development, and this yields a 
relation between (A, to, h): 



TV 



1 

dx— 

-i (A 



VT 



2x) 2 + u 2 



(25) 
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Figure 2: The transition line h*(uj), estimated by analytical, numerical and heuristic 
tools, together with the small frequency behavior h* oc u; - 25 . 



from which we estimate the transition line where the time scale t r (u,h) given in 
Eq.(|21|) diverges. This is equivalent to the condition \(u, h*) = 2. It is easily 
shown that in the small frequency limit, the scaling t r °(h*) ~ r is recovered, in 
accordance with our heuristic argument, Eq. fl2"2"|). In Fig. ^|, the analytic estimate 
for the transition line h*(u) is compared with its direct numerical estimation, and 
with our heuristic argument. Also plotted is the small frequency behavior of the 
critical field, h* oc u - 25 . 

Moreover, Equation ( ppp implies 

lim A(u;, h) = X 

which means that, in the limit uo — > 0, the field acts as a constant field with an 
amplitude given by its root mean square value, which is physically reasonable. Im- 
portantly enough, this means that the effect of a dc field with uj strictly zero and 
the limit of an ac field with vanishing frequency are different. In the case p = 2 this 
feature has no effect in the value of the transition field h* since it vanishes in both 
cases. However, when p > 3 we shall find a non trivial consequence of this result 
(see Section f|). 

3.3 Below the transition h < h*(u): aging in an ac field. 

Below the transition, Eq. ( ^0|) shows that the energy density slowly converges to- 
wards its asymptotic value e^ui, h < h*) = — 1 as a power law e(t) — ~ t^ 1 with 
a superimposed oscillation at a frequency u/n. This asymptotic value is indepen- 
dent of h. A slow (here a power law) convergence of one time quantities is typical 
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Figure 3: Correlation function of the p = 2 model atT = 0for/i = 0.5,c<j = l and 
waiting times t w = 5, 10, 20, from bottom to top. The points fall on the discrete 
times t = nr, with n an integer. 



of aging systems [11]. The angular velocity 2oj is expected since Eq. (^J) remains 
unchanged by the transformation t ^ t + tt/uj; — > — s M ; z — > z. The constraint 
z(t), and the energy density, are then n/u-per iodic. 

That tti is also lost is demonstrated by looking at two-time quantities, typically 
correlation functions. The behavior of C(t,t w ) is represented in Fig. || This figure 
also shows the interesting feature that the dynamics can be decomposed into two 
well-separated time scales. At short time separation t — t w , the time-scale for the 
approach to the plateau at C — q EA (note that q EA = 1 at T = 0) does not depend 
on t w , whereas the decay from the plateau towards zero arises in a second time scale 
that clearly depends on t w . During the waiting-time dependent decay the curves 
have oscillations. Quantitatively, the behavior of the correlation function is entirely 
dominated at long times by the term C(t,t w ) ~ f (^ ± ^ L ) / (il(t)n(t w )) . This implies 
that, apart from the oscillation, the correlation scales as t/t w for t ^> t w ] we have 
numerically checked this point. 

In tapping experiments, data are obtained for times of the form t = nr, where n is 
an integer. We have represented these data by points in Fig. || and the experimental 
measurements would be very similar to the usual aging case. 

The response function is analyzed in Fig. [|, where the parametric plot of the 
integrated response against the correlation is built. As expected, there are strong 
violations of the fdt. These violations are however similar to those encountered in 
the relaxational case, in the sense that there are two distinct behaviors, depending 
on the time scale considered. For a time separation t — t w <^t w , fdt is not satisfied, 
and no effective temperature can be defined. For the second time scale t > t w and 
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Figure 4: Parametric plot of the integrated response-correlation of the p = 2 model 
at T = for h = 0.5, u = 1 and waiting times t w = 5, 10 and 20. For the fdt to be 
satisfied, the points have to lie on the vertical line. The inset represents the same 
data, but for times t = nr; the symbols are the same as in Fig. 



longer, on the contrary, all points fall nicely on a straight line, with a superimposed 
oscillation. The bound (Hp shows however that the deviation from FDT can only 
be significant at large times, when both dC(t,t w )/dt w and R(t,t w ), integrated over 
a period, are small. This is confirmed by Fig. |j. 

The inset of Fig. ^ shows x(C) for the discrete times t = nr. If one uses 
this stroboscopic measurement, it remains clear that FDT modifications are well 
accounted for by the ansatz fllpp . Analytically, neglecting as above the second 
term in the correlation function, the effective temperature defined through Eq. (|9|) 
scales as T EFF ~ tw — > oo. An infinite effective temperature is also present in the 
relaxational dynamics of this model, and is typical of domain growth models. 

With these simple considerations we have argued that the effect of a small ac field 
does not change the aging behavior of this model at T = 0. This is to be confronted 
to the effect of non-Hamiltonian perturbations that change the aging scaling, as 
seen for instance in the correlation function that becomes C(t,t w ) ~ f(exp(y/t^ — 
v£)) HI8| . At non-zero temperature, shear-like forces introduce a finite relaxation 
time after which the dynamics is stationary. We have not solved analytically the 
p = 2 model at finite temperature in an ac field, but the numerical solution of the 
full equations suggests that there is a finite region of the phase diagram (T, h, uj) 
in which aging effects survive. This explicitly shows that the effect of shear and 
tapping are rather different. 
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3.4 Above the transition: driven steady state. 

Above the transition, the situation is simpler. After a time scale t r (u, h) = [X(uj,h) — 
2] _1 the slow dynamics is lost. The relation ( p5|) implicitly determines A, which is an 
increasing function of the amplitude h of the driving force. In the stationary state, 
the energy density oscillates with an angular velocity 2u around its asymptotic value 
eoo(u),h > h*) = — \(u>, h)/2. The behavior of the correlation is dominated by the 
second term in Eq. fllTf). This shows that it decays towards zero in a time scale of 
order t r (u, h): there is no more aging. 



4 Behavior for p > 3 



Going beyond the solution of the unperturbed case |5(J to solve the set of cou- 



pled integro-differential equations dJ) analytically is a very hard task. In order to 
illustrate the main properties of the solution, we solved Eqs. (^) numerically, by 
constructing the two-time solution step by step in time. We first demonstrate that 
also in the case p > 3, there exists a transition line below which glassy properties of 
the undriven model persist, and study then the fdt- violations in both phases. 

4.1 Evidences for a dynamical transition 



Let us recall first the effect of a dc magnetic field of amplitude h on the system |63 
In contrast to the spherical Sherrington-Kirkpatrick, the spin glass phase may exist 
in a dc magnetic field. There are both a dynamic and a static transition between 
the spin glass and the paramagnetic phases. At T = 0.2, the dynamic transition 
takes place at h*(u> — 0) ~ 1.1 [B3 



As in the case p = 2, we begin our study by focusing on the behavior of the energy 
density e{t). It is displayed in Fig. || at temperature T = 0.6 and angular velocity 
<jj — 1, for different field amplitudes. The long time behavior of e(t) indicates that 
the transition occurs for a field 0.6 < h* < 1.0 at this frequency. The asymptotic 
value of the energy density slightly increases below the transition, and reaches a 
much larger value, when h > h*. 

In the following we concentrate on the temperature T = 0.2 ~ 0.32 T c , i.e. well 
below the zero-field transition. We focus first on the dependence of the two-time 
autocorrelation functions on the amplitude of the applied field. In Fig ^|-a, the field 
is h = 0, and we observe the usual aging f5"0 |. It is then clear that aging is still 



present if the field is not too strong, h — 0.1 and h — 1 in Figs. ^)-b and c, whereas 
at a stronger field, h = 2 in Fig. |6|-d, the correlation very rapidly tends to zero, with 
a superimposed oscillation. These observations reinforce the evidence in favour of 
a dynamic transition which, at T = 0.2, occurs at a field strength 1 < h* < 2, for 

LU = l. 

We turn now to the question of identifying the critical line in the u direction. 
Fig. [7| shows the auto correlation function at fixed amplitude of the applied field and 
several values of the angular velocity. In Fig. |7|-a, the frequency is zero, and aging 
is observed with a two-step decay of the correlation. Note that the second decay 
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Figure 5: Time dependence of the energy density under an ac field with amplitude 
h = 0,0.3,0.6, 1 (from bottom to top) and angular velocity to — 1. The curve for 
h — 1 approaches rapidly an asymptote, whereas e(t) still evolves at large times for 
the other cases. 



is towards a value q Q > 0, as opposed to the zero-field case f63j. The frequency in 
increased in Figs. 0-b,c,d. At intermediate frequency, Fig. 0-b, aging is suppressed. 
When ui is further increased, Figs. 0-c,d, it is clear that the glassy phase is reentrant. 

It is of course very difficult to determine the transition line h*(u) with accuracy, 
by the numerical solution of the dynamical equations. We are hence not able to draw 
a figure similar to Fig. ^| for p > 3. We show in Fig. |^ a schematic representation 
of this critical line with the feature, already encountered for p = 2, that the limit 
uj —>■ is peculiar. From the numerical solution at uj > it seems that 

lim h*(T, uj) < h*(T, a = 0) . (27) 

This is represented in Fig. |8] and explains the reentrance of the glassy phase described 
in Fig. 0, when the frequency increases. Note however that we cannot numerically 
discard the possibility of a non-monotonous, but smooth, behavior of the critical 
line at u < 0.1. 

4.2 Fluctuation-dissipation theorem 

We turn now to the study of the fdt. In Fig. [5] we show the evolution of the x vs C 
plots with t w . In Fig. |9|-a, the system is not in the glassy phase. Strong violations 
of the fdt are observed, but no effective temperature can be defined. The system 
is athermal. 

In Fig. ^|-b, on the contrary, the system is in its glassy phase. One recovers then 
a x{C) curve that is very similar to the one in zero field. The first part is almost a 
straight line with slope — 1/T while the second decay follows a temperature — 1/T EFF . 
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Figure 6: Each panel shows the auto-correlation C(t + t w ,t w ) as a function of t, for 
four values of the waiting-time t w = 3, 12,49, 198 at T = 0.2. The angular velocity 
of the applied field is u = 1. In panels a to c the scale is logarithmic, the amplitudes 
of the applied fields are h = 0, 0.1, 1 and the system is in its glassy phase. In panel 
d the scale is linear, h = 2 and the slow dynamics is suppressed. 

As mentioned in Section 3, the fact that for small time scales, the fdt is nearly 
satisfied results directly from the bound ([12]). 

The non-trivial outcome of this study is the fact that a well-defined T EFF , which 
exists in the non-driven system, may still be defined below the dynamical transition, 
in presence of the drive. A stroboscopic construction, as the one in Section 3.3, will 
yield a perfect straight line that defines T EFF unambiguously. 

The actual value of the effective temperature may depend on various parameters. 
We have searched for dependences on the frequency and amplitude of the ac field. 
Within the range of parameters we could numerically explore, we have not observed 
large dependences. Figure [l^ display the \ vs C curves for two choices of field 
amplitudes, and for two choices of angular velocity at fixed field amplitude: the 
dependence is indeed weak. 

4.3 Beyond mean-field: Effect of trapping states 

In this Section we present results from a numerical simulation of the Ising version 
of the p-spin model in an ac field, in the particular case p = 3. We focus on 
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Figure 7: Each panel shows the auto-correlation C(t + t w , t w ) as a function of t, for 
four values of the waiting-time t w = 3, 12,49, 198 at T = 0.2. The strength of the 
applied field is h — 1 and the frequencies are uj = 0,0.1,1,10 from panel a to panel 
d. The transition in an ac-field of amplitude h — 1, from the liquid to the glassy 
phase, occurs at a frequency 0.1 < u c < 1. 




Figure 8: The full line is a schematic representation of the transition line h*(u) for 
p > 3. The angular velocity in the four panels in Fig. [7|move from left to right along 
the horizontal dashed line (see the text for more details). 
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Figure 9: The x(C) curves at T = 0.2 with an applied ac-field of strength h = 1 and 
frequencies to = 0.1 (liquid phase) and to — 10 (glassy phase), for different waiting 
times t w = 6,12,24,49. 
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Figure 10: The x(C) curves at T = 0.2. Panel a, the two curves correspond to h = 
and h = 1, uj = 1. In both cases the system is in its glassy phase. In panel b the 
two curves correspond to uj = 1 and uj = 10, both with an amplitude h = 1, The 
lines are guides to the eye and they have the slope of the second part of the decay. 
t w = 24 for all curves displayed. 



the time dependence of the energy density e{t) = N^ 1 J2i<j<k JijkSiSjSk and the 
magnetization density m(t) = A^~ 1 ^ i Sj(t). The sizes N = 50 and 150 have been 
used, together with a constant temperature T = 0.01, in all simulations. The 
amplitude and the frequency of the magnetic field have been varied. 

The interest of such an investigation is that the finite N behavior of the model 
is accessible. The system is hence able to escape and visit the trapping states 
described in Section |2] in a finite times |Tj|. For that reason, the results presented 



in this section cannot be obtained in a mean-field approach. 

In the absence of magnetic field, the energy density rapidly approaches its asymp- 
totic value. At such a low temperature, the system reaches a metastable stable, 
characterized by a nearly zero magnetization density, and cannot escape within the 
numerical time window. Fig. [TT] (left) illustrates this 'jammed' behavior. The value 
of the energy density depends on the initial conditions, which means that the sys- 
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Figure 11: Evolution of the energy (bottom curve) and magnetization (top curve) 
densities without (left) and with (right) magnetic field. The parameters of the field 
are h — 2, to — 0.01. Note that the time range in the left figure is much shorter, 
since both e(t) and m(t) are constant for times t > 50. 



tern may be blocked in different trapping states when it starts from different initial 
conditions. 

The influence of an ac magnetic field on this jammed behavior is evident in 



Fig. [IT]. Both the energy and the magnetization density are oscillating functions, 
and their behavior is intermittent. For some time windows, both quantities have an 
evolution which varies considerably from period to period. For other time windows, 
on the contrary, the different cycles are very similar. This can be interpreted as 
being due to the presence of the trapping ('jammed') states: the system usually 
evolves inside one of the trapping states, without escaping, and all the periods are 
equivalent. But from time to time, the driving force is able to make the system 
escape the state, and the evolution is very erratic until another trapping state is 
found. 



This is confirmed in Fig. [12|, where the energy density is represented as a function 
of t for times of the form t = nr. The long horizontal plateaux are the moments 
where the system is jammed, whereas between the plateaux, the energy density 

we show in Fig. [1^ that the time 
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changes values very rapidly. As in Ref. 
evolution of e(t) is 'self-similar', in the sense that zooming on a time window makes 
smaller plateaux become visible while the overall evolution looks the same. 

This intermittent behavior is clearly reminiscent of the behavior encountered in 
some granular experiments, where the powder is very slowly perturbed, like the ones 
performed by the Jussieu group [Q. In particular, it would be very interesting to 
perform a more careful analysis of the statistical properties of e(t) like, e.g., mea- 
suring the statistics of trapping times p3 |. We note also that this phenomenology 
is very similar to the one of the so-called 'trap model' []65|1 , that has been extended 



by Head to describe the phenomenology of granular materials 



Finally, it is interesting to stress that even within one cycle of the field, the 
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Figure 12: Same parameters as in Fig. |IT]. Evolution of the energy density with a 
magnetic field. Only one point per period is represented. 



evolution is not regular at all. This can be nicely seen in Fig. [13], where the field 
h(t) is represented as a function of the magnetization m(t). This is the usual view of 
a hysteresis loop in ferromagnetic systems. It is clear that the shape of the loops is 
far from elliptic, and that the system evolves in steps, rather than continuously: this 
is analogous to the Barkhausen noise [EfJ. The overall shape of a cycle drastically 
depends on the fact that the system is trapped or not. It is beyond the scope of this 
paper to study these loops in detail, but we emphasize that this (mean-field) model 
could be an interesting starting point to study the Barkhausen noise, in the spirit 
of Ref. 



5 Summary and link with granular matter 

Phase diagram. - Our findings concerning the behavior of glassy systems under 
a time-dependent driving force are summarized in Figs. |8| and [14]. For any fixed 
frequency of the applied field, and a temperature T < T c , where T c is the zero-drive 
glass transition temperature, there are two regimes: (I) At small drive, the system 
exhibits slow ("glassy") dynamics. (II) At large drive, there is no slow dynamics. 
There exists then a well-defined critical drive separating these two regimes. On 
the other hand, if we work at fixed drive and modify the angular velocity, the 
system undergoes a transition from a liquid-like phase at low (though non vanishing) 
frequency to a glass-like phase at high frequency. This result is displayed in Fig. || 
The phase diagram in Fig. |TJ| is the analog to the (T, T) plane of the three- 
dimensional phase diagram proposed by Liu and Nagel 0. We have then inves- 
tigated its properties in detail, which had not been done so far in the case of a 
time-dependent driving force. Contrary to previous phenomenological modeling of 
the dynamics of granular matter assuming that the drive T is related (in a possibly 
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Figure 13: Hysteresis loops. Panels a and b correspond to the period of rapid motion 
while panels c and d correspond to the interval when the system is trapped. N = 50, 
T = 0.01, h = 2 and u = 0.1 



nonlinear way) to the temperature T of a glassy model borrowed from statistical 
mechanics, we propose here to study precisely the interplay between T and T as an 
intermediate step from glassy to granular materials. The conclusions we draw are 
then directly pertinent to gently driven granular materials or slow granular rheology. 

We emphasize once again, that the model we have studied is not intended to 
describe granular matter in full microscopic detail, but rather to act as a source 
of inspiration for the interpretation of numerical and experimental results. It may 
also motivate new experimental measurements. Let us however briefly summarize 
its dynamic behavior in the context of existing data for granular matter. 

One-time quantities. - In phase (I), one-time quantities (we have focused on 
the energy density) typically decay towards their asymptotic values as power laws. 

It is experimentally established that the density p(t) in a gently driven granular 
system exhibits a very slow relaxation IJ. The parallel can be drawn with the 
energy density e(t) (or more precisely with — l/e(t)) in glassy models. The seminal 
experiments at Chicago have exhibited a logarithmic relaxation of the density @], 
as found in all subsequent works J|, ^6], ^7], |28|, ^9], 
Ref. [§] explicitly excludes a power law behavior.} 
of Nicolas et al. [EJ showed that the logarithmic law is not a universal behavior 



§11], 0, [34], |35|, [3§]. (Note that 
However, the recent experiment 
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Figure 14: Temperature-Drive phase diagram of a glassy system. The vertical axis 
represents the typical amplitude of the driving force: its precise definition depends on 
the specific experiment considered; T c is the zero-drive glass transition temperature. 



and that the time dependence may depend on the way the system is perturbed. 
From the theoretical point of view, one could prefer a power law decay for any 
one-time quantity since it does not involve an intrinsic time scale, whereas the 
logarithmic law l/ln(t/to) explicitly introduces a time scale to- in this context, it is 
interesting to note that Barrat and Loreto [p8|| , in their study of the Tetris model, 
have shown that the fitting parameter t actually depends on the considered time 



window; qualitatively t oc t w . Head [36] had already emphasized the role of the 



time window on the fitting parameters, in particular on the asymptotic value of the 
density. 

In the context of glassy systems, it is a well-known feature that the interpretation 
of data over a finite time-window is always delicate. A typical example is provided 
by the case of dipolar glasses where the dielectric constant for samples prepared 
with different cooling rates seems to approach different asymptotic values [|1^ . 



Two-time quantities. - In phase (I), the model we have studied also exhibits 
aging. For p = 2, we find a t/t^-scaling (in stroboscopic time) of the autocorrela- 
tion function while for p > 3 a much more careful numerical analysis is needed to 
determine this law. Other glassy models perturbed with ac-forces may lead to other 
aging scalings. Generally speaking, theoretical arguments |TT| indicate that one can 
expect the following behavior for any two-time function F(t,t w ) 

F(t,t w )~r(^L), (28) 



h{t, 

in a given time-scale. T and h are two model/system-dependent scaling functions. 
More complicated scalings, such as a sum of different terms like Eq. (p8|) implying 
several different time scales in the problem, are also possible. 

Two-time quantities such as the density-density correlation function have been 
numerically studied, clearly demonstrating that tti is lost P7L |3~I], EH|: the Tetris 
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model exhibits a ln(t) / \n(t w ) behavior [P? , p8f , whereas the parking lot model shows 
at/t w aging behavior [ plfl . 

An experimental determination of the scaling of two-time quantities will be use- 
ful to discriminate among the different models proposed to describe granular matter. 



This has already been proposed in Refs. [27, EM Rl| and amounts to an experimental 



determination of the function h(t) involved in Eq. (p8|). Such experimental deter- 
minations in glassy systems are by now numerous. They have emphasized on the 
one hand the universality of the aging phenomena in many microscopically different 



systems [[12], [13|, [14|, |15|, pT6[| , but on the other hand they have also revealed subtle 
differences between them [|TT], |T3"| . 

Fluctuation- dissipation theorem violations. — Our study has shown that modi- 
fications of the fdt, rather similar to those observed for aging glassy systems, also 
arise in glassy models driven by time-dependent forces. The parallel with granular 
materials suggests that this may also happen in driven powders. 

The parametric plot we have built clearly shows that driven glassy systems have 
two distinct time scales, which we interpret as in Refs. || |§. A fast one, which 
is essentially independent of the waiting time, represents the individual motion of 
grains: no effective temperature, in the sense of Eq. @, can be here defined. Since 
the model we have used is itself thermal (it uses a Langevin equation with a thermal 
noise) the fast motion has a reminiscence of the temperature of the bath (that 
may be zero). In a more realistic model of granular matter, it is likely that the 
fast motions will be completely athermal. On the contrary, the long time scale 
growing with the waiting time and giving rise to the t/t^-relaxation, represents 
large structural rearrangements in the powder. Our results indicate that these slow 
degrees of freedom are thermalized at a well-defined effective temperature. We are 
confident that this result also holds in realistic models of granular matter. 

Again, numerical and experimental work could be used to check these predictions 
in granular materials. We have shown that the two models studied here have differ- 
ent fdt violations: the p = 2 case has an infinite effective temperature, whereas the 
p > 3 models have a single finite effective temperature. Other glassy models (for 
instance mean-field spin glasses) have been shown to have a hierarchy of effective 
temperatures when they are undriven, and shall probably keep this behavior under 



driving forces |h||. We conjecture that granular materials fall in the p > 3 category, 
by analogy with structural glasses. This prediction has to be experimentally and 
numerically tested. It is not unlikely that different microscopic models lead to dif- 
ferent FDT-violations; such a numerical determination could once again discriminate 
between them. 

To study FDT violations during the compaction experiment, one has to com- 
pute, separately, the correlation function between two observables and its conjugated 
response function. In experiments, for instance, the Nyquist relation is checked 



through dielectric measurements |14|, |68[. In granular materials, a natural choice is 
to check the Einstein relation between self-diffusion and mobility of the grains. The 
results obtained in Refs. |69|, |70| are then very encouraging, although the driv- 
ing field is not time-dependent in these studies. Nicodemi studied instead the 
fluctuation-dissipation relation between height-height displacement and its conju- 
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gated integrated response to small shaking amplitude perturbations in the Tetris 
model fTIJ . In the low density regime (fluide-like), after a short transient, these two 
quantities are related by the usual FDT. In the high density regime, the relation 
between these quantities shows much stronger deviations from FDT than the model 
we studied here. This difference may be due to the kind of perturbation used by 
Nicodemi. For a thermal system, the procedure in Ref. [FT] is indeed similar to 



measuring the energy response to a temperature change: The resulting FDT prop- 
erties are found to be of an unusual type [72] also in this case. However, the notion 
of an effective temperature [52j implies that the FDT relations of all observables 
evolving in the same time-scale should be identical. As suggested in Refs. [2R IBS 



an interesting open question is: does a measurement in the Tetris model done with 
a more standard infinitesimal perturbation lead to the same FDT relation? This 
point deserves a more detailed investigation in this and other glassy models. 

The dynamic transition, and the regime (II). — We have found a well-defined 
transition in the phase diagram of Fig. [TJ]. This is in complete agreement with ex- 
perimental and numerical observations of the existence of a critical value of the drive 
r*, below which glassy effects may be observed. The existence of this dynamical 
transition provides a very natural context for the interpretation of the numerous 
non-equilibrium effects encountered in granular matter. Indeed, the transition line 
in the Fig. |14] links the standard glass transition arising at T = T c , T = in the 
glassy model to the dynamical transition arising at T = 0, T = T* in the driven 
dynamics of athermal systems. In our opinion, this result gives a nice theoretical 
support to the glasses/granular analogy, making a deep connection between the two 
situations. 

We wish to emphasize that it was not evident a priori that the transition could 
survive a finite time-dependent driving force. This is one of the main results of 
this paper and it has to be confronted to the very different effect of shear-like forces. 
Indeed, an infinitesimal shear-like perturbation, is enough to introduce a finite time- 



scale in this model, and aging is hence interrupted |T8l 21 



The liquid regime (II) is less interesting in the present context, since no slow 
dynamics is present. Let us stress that the system is still strongly driven, and hence 
no 'equilibrium' state, in the sense of statistical mechanics, is reached. In particular, 
no effective FDT temperature can be defined: the system is completely athermal. 
The very existence of an effective temperature relies indeed on the existence of two 
well-separated time scales allowing the slow degrees of freedom to thermalize and 
the susceptibility / correlation parametric curves of Fig. ^|-a are not very useful to 
describe this non-equilibrium situation. 



6 Conclusion: towards a thermodynamical descrip- 
tion of the slow granular rheology? 

The interpretation of FDT violations in terms of an effective temperature T EFF has 
been developed in Ref. ||52|| . Having a well-defined concept of temperature is a 
crucial first step for a thermodynamical description of granular materials [[7J3] . It is 
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then important to discuss the possibility of a link with previous thermodynamical 
concepts in the granular literature. 

The notion of a 'granular temperature' has been introduced to extend thermody- 
namics to strongly perturbed powders [|l| . In analogy with the kinetic gas theory, it 
has been assumed that distribution of the grain velocity v is Maxwellian and via the 
equipartition relation T oc (v 2 ). We have been concerned with the weakly perturbed 
regime, where the dynamics is slow. This hydrodynamic definition is not supposed 
to be relevant in this case. This is also well known in the field of glasses where 
nontrivial effective temperatures are known to exist while the kinetic energy is an 
observable that very rapidly equilibrates with its environment, leading to a kinetic 
temperature that coincides with the temperature of the bath. Very recently ]74j] , an 
attempt to extend the hydrodynamic theory to a regime where the grains are not 
fluidized has been proposed, also using the concept of kinetic granular temperature. 
It would then be very interesting to try to understand the two concepts in a unified 
way. 

More related to ours is the approach of Mehta et al. |§, discussed in the in- 
troduction. This phenomenological two-step, two-temperature model finds a nice 
justification within the scenario emerging from our results. This model was in fact 
inspired by the illuminating work initiated in the late 80's by Edwards and cowork- 
|38[| . They have postulated thermodynamic relations in analogy with the usual 
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'thermal' thermodynamics, where the volume V plays the role of the internal energy 
U. In this approach, the entropy is the logarithm of the number of configurations 
with volume V . The so-called compactivity, X, is the analog of the temperature 



and it then defined as 5c 



- = — (29) 
X ~ dV ' 1 ] 

still by analogy with the usual definition T~ l = OS/ dU . The connection between the 
compactivity X and the fdt temperature T EFF discussed in this paper has already 
been explored ]10|, [57], [58], |62|, [7I| . Crucially for our study, Edwards 7 definition of 



the compactivity coincides with the (asymptotic) fdt effective temperature defined 



in Eq. (14). This result holds for mean-field models like the one discussed here [62 



as well as in the lattice gas model studied in Ref. ^7j. However, these studies have 



been done without an external forcing. It would be very interesting to extend the 



numerical analysis of Ref. |57j to the same model now driven by an ac force. 



We have obtained here that the existence of an effective temperature for the slow 



decay resists a finite ac force. This supports the conjecture that the definitions ( |T4| ) 
and (^) may be of fundamental interest for the study of glassy /granular materials. 
Computing them in realistic models is then a challenge for future research [^, ^ 



In this respect, very recent works |76j studying spin models on random 



graphs may give some insights on the role played by metastable states, which can 
have an important influence on the dynamics, as we showed in Section [O 



In conclusion, we have found that the study of driven glassy systems provides a 
theoretical framework to understand the slow granular rheology. The existence of 
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a dynamical transition justifies well the use of 'modified glassy models' to describe 
granular materials. The existence of an effective temperature for the slow degrees of 
freedom provides in particular a nice theoretical basis to previous 'thermal' models 
for granular matter, and to the seminal approach of Edwards. 
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